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Abstract 

We examine the role of ^ (T7r° ^ 7r’'"7r“7r'^ decay in the Dalitz plot analysis of ^ pn ^ 

7r^7r“7r° decays, employed to extract the CKM parameter a. The (ttf channel is significant because it can 
break the relationship between the pengnin contributions in —> p^if', B —> p''"7r“, and B p“7r+ decays 
consequent to an assumption of isospin symmetry. Its presence thns mimics the effect of isospin violation. 
The state is of definite CP, however; we demonstrate that the B ^ p-K analysis can be generalized to 
include this channel without difficulty. The a or /o(400 — 1200) “meson” is a broad 1 = J = 0 enhancement 
driven by strong tttt rescattering; a suitable scalar form factor is constrained by the chiral dynamics of low- 
energy hadron-hadron interactions — it is rather different from the relativistic Breit-Wigner form adopted in 
earlier B —> <jtt and D on analyses. We show that the use of this scalar form factor leads to an improved 
theoretical understanding of the measured ratio Br(i3'’ —> p^7r*)/Br(i3~ —> p°n~). 
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1 Introduction 


Measurements at SLAC and KEK of the time-dependent CP-violating asymmetry in B{B) —> J/ipKs [Q, ||, 
yielding sin(2/3), have conclusively established the existence of CP violation in the B meson system. The results 
found are consistent with Standard Model (SM) expectations Q, so that establishing whether the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix Q is the only source of CP violation in nature, as in the SM, or not requires 
the empirical measurement of all the angles of the unitarity triangle. 

In this paper we consider the determination of a through a Dalitz plot analysis of the decays B^{B^) 
piT 7r+7r“7r° under the assumption of isospin symmetry a I- Ten parameters appear in the analysis, and 
they can be determined in a fit to the data. Nevertheless, the assumption of p dominance in B Sir decays 
has no strong theoretical basis Q, so that the contributions from other resonances in the pir phase space may 
be important. We discuss how the isospin analysis can be enlarged to include the cttt channel as well. The 
cr or /o(400 — 1200) “meson” is a broad J = I = 0 enhancement, close to the p meson in mass, so that the 
CTTT channel can potentially populate the Stt phase space associated with the pir channels. The (Ttf final state 
contributes preferentially to the p°7r° final state. In the context of the isospin analysis, such contributions are 
of consequence as they invalidate the underlying assumptions of the isospin analysis and thus mimic the effect 
of isospin violation. 

Our considerations are inspired in part by recent studies of D~ —> 7r“7r+7r“ decay: the E791 collaboration 
find that the pathway D~ —> —*■ 7r“7r+7r“ accounts for approximately half of all D~ —> 7r“7r+7r“ decays 

Deandrea and Polosa have argued as a consequence that the B —> cttt channel contributes significantly to the pir 
phase space in i? ^ tttt+tt” and modifies the ratio B{B^ —> p'^Tr^)/B{B~ p°7r“) to yield better agreement 

with experiment |^. The scalar form factor, which describes the appearance of the a in the tt+tt” final state, 
enters as a crucial ingredient in the assessment of the size of these effects. The scalar form factor cannot be 
determined directly from experiment; nevertheless, ample indirect constraints exist, permitting us to describe 
its features with confidence |T^ . Nevertheless, different approaches, with different dynamical assumptions, yield 
roughly comparable descriptions of the tttt scattering data, so that the emergence of a favored form of the scalar 
form factor does not resolve the question of whether the cr is a pre-existing resonance or, rather, a dynamical 
consequence of tttt interactions in the final-state. We follow Ref. [Q and adopt an unitarized, coupled-channel 
approach to the final-state interactions (ESI) in the tttt-KK system, and match the resulting scalar form factor 
to chiral perturbation theory (CHPT) in the regime where the latter is applicable. The resulting form factor, 
in the tttt channel, discussed in Sec. 5, is strikingly different from the relativistic Breit-Wigner form adopted 
by the E791 collaboration in their analysis of the tr in £>+ ^ tt+tt+tt" decay — the latter form factor is also 
used in Ref. Q. The differences are particularly large as ^/s —> so that the relativistic Breit-Wigner form 

is at odds with CHPT in the precise region where it is applicable, note Fig. 4 in Ref. 0. This casts doubt on 
the recent conclusions of Refs. |^, |^, prompting new analyses incorporating a suitable scalar form factor. 

The generation of the a resonance via strong rescattering effects, as in the approach we adopt, indicates 
that OZI-violating effects in the scalar sector are significant. Moreover, the “doubly” OZI-violating form factor 
(0|ss|7r7r) is non-trivial as well; such a contribution is needed to fit the tttt and KK invariant mass distributions 
in J/ip (j)7nT{KK) decay |]^. These effects are also needed to explain the branching ratios of the decays 
of the ao(980) and /o(980) states into tttt and KK final states [0. These observations give new insight on 
rescattering effects in hadronic B decays, generating a new mechanism of factorization breaking in n > 3 particle 
final states. 

The contribution of the B —> cttf channel to the B —> p^i: phase space can also modify the inferred empirical 
branching ratios in these channels. Combining the CLEO results 

Br(R-^p°7r-) = (lO.dtfj ± 2.1) • 10"® , (1) 

Br(R° ^ p±7r=F) = (27.6± 4.2) • 10"® (2) 

with the BaBar result |0 Br(i?° ^ = 28.9 ± 5.4 ± 4.3 (charge conjugate modes are implied) yields, 

adding the errors in quadrature and ignoring correlations. 


Br(R° ^ 

Br(R“ ^ p^'K~) 


= 2.7± 1.2 . 


( 3 ) 


This ratio of ratios is roughly 6 if one works at tree level and uses the naive factorization approximation for 
the hadronic matrix elements |p^. The inclusion of penguin contributions can alter this result, and potentially 
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yield better accord with theory and experiment ||^, |^, |^. However, our focus will parallel that of Ref. [§: 

we wish to examine how B —> (ttf —> Stt decay, given a particular scalar form factor, can effectively modify the 
theoretical prediction of the ratio given in Eq. (i- It is apparent that B —> cttt is of greater impact in H —> p°7r 
decay, so that the inclusion of such contributions ought alter the ratio of ratios. 

We begin by reviewing the isospin analysis in B^{B^) ^ p7r —> 7r+7r“7r° decay [^ ||] in Sec. ||, and discuss its 
extension to include B°{B°) —*■ cttf ^ 7r+7r“7r°) decay in Sec. ||. We proceed by evaluating cr-mediated i? —> Stt 
decay in Sec. |[ relegating the p-mediated i? ^ Stt decay formulae to App.^ Our analysis employs the scalar 
and vector form factors discussed in Sec. || and Sec. respectively. We conclude with a presentation of our 
results in Sec. and an accompanying summary. 


2 Preliminaries: Isospin analysis oi B —^ pn 

Let us recall the isospin analysis possible in H ^ p7r decay il- Under the assumption of isospin symmetry, a 
pTT final state can have isospin If = 0,1, or 2, whereas the i?+, B^ states form an isospin doublet. Thus we can 
have |A/| = 1/2, 3/2, or 5/2 transitions in i? —> p7r decay, so that we can parametrize the amplitudes which 
appear by A|A 7 |, 7 j. We have |^, 


a+_ = A(R° 


P+TT ) = 


1 


a_+ = A(R° 


■) = 


and 


Qoo = 


2V3 

1 

2^3 

„o„o 


[^3/ 


2,2 


■A 


5 / 2 , 2 ] 


[^3/2,: 


■A 


5/2,2l 


[^3/2,1 

[^3/2,1 


■ ^1/2 
■ Ai/2 




5 " ) - --^[^ 3 / 2,2 


U/2,2i 


^^/2.0 


( 4 ) 
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( 6 ) 


noting that A(B^ 7r+7r“7r°) = f+a-i _+ /_ a_+ /oaoo, where fi is the form factor describing p* ^ 

TTTT. Isospin is merely an approximate symmetry of the SM; nevertheless, our parametrization possesses three 
independent “isospin” amplitudes, distinguished by If, to describe the three empirical amplitudes aij, so that 
it persists in the presence of isospin breaking as well. 

The lowest-dimension operators of the effective, |Ai?| = 1 Hamiltonian generate transitions of |A/| = 1/2 
or 3/2 character, so that a |AJ| =5/2 transition is generated, in this order, through long-distance, isospin¬ 
breaking effects in concert with a |A/| = 1/2 or 3/2, short-distance, weak transition. If we neglect transitions of 
|AJ| =5/2 character and, indeed, isospin-violating effects all together, the transition b —!■ qqd which mediates 
B ^ pn decay can be realized through a “tree” amplitude with |A/| = 1/2 or 3/2 or through a “penguin” 
amplitude with |A/| = 1/2. Practically, the decay topologies are distinguished by their weak phase, so that the 
contributions associated with the CKM factors U^bl^d’ defined to be tree contributions, regardless of 

their dynamical origin. The unitarity of the CKM matrix in the SM implies that two combinations of CKM 
factors suffice in describing b —>■ qqq'] here we associate the combination with the penguin contribution. 

Noting 


v:bVud 

KbVudl 


= e" 


VttVM 

KbVJ 


= 6“*^ 


( 7 ) 


and a = TT — /3 — 7 , we have 


e*^a+_ 

e*^a_^ 

e*^aoo 


= T+- 


+ P+- 
+ P-+ 


T-+e““ 


rpOO ^ — _j_ ^00 


( 8 ) 


The overall weak phase is without physical impact and can be neglected, because in the SM the weak phase 
associated with mixing is controlled by q/p = exp(— 21/3)^. Thus qdij/p cx exp(—i/3), just as Oy 

is. Consequently, the isospin analysis in B ^ pn decay determines a. The crucial assumption of the isospin 

flip the overall sign of the apo amplitude with respect to Ref. ||Ie|], to conform with our computation of the amplitudes. 

^"^Recall that the B mass eigenstates are defined via \Bij) = p\B^) + q\B^) and \Bh) = p\B^) — q\B^). We assume throughout 
that the width difference of the two B mass eigenstates is negligible, so that \q/p\ = 1. 
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analysis is to associate the CKM factor with |A/| = 1/2 transitions exclusively, so that from Eqs. 

we have 

pOO = i(P+- +p-+) . (9) 

The overall strong phase in Eq. is trivial, so that with Eq. we have ten parameters in all, which can be 
determined in an analysis of the Dalitz plot |^, . 

The presence of the (ttf final state in the phase space associated with the p°7r° channel breaks the relation 
assumed in Eq. (|^), and thus mimics the appearance of isospin violation. In this paper we study how the impact 
of this additional decay channel can be minimized. It is worth noting, however, that the utt final state is of 
definite CP, so that the isospin analysis can be enlarged to include this channel as well — additional observables 
are also present in this case. Before doing this, let us enumerate the ways in which SM isospin violation can 
impact the usual B pir analysis, to determine whether the impact of these effects can be reduced as well: 

i) Isospin violation can generate an additional amplitude, of |AJ| = 5/2 character, as in Eqs. A 

|A/| = 5/2 amplitude can be generated by 0{md — ruu) or 0(a) effects in concert with a |A/| = 3/2 
weak transition, or by 0{a) effects in concert with a |A/| = 1/2 weak transition. The 0{md — rriu) term 
acts as an isovector interaction. We recall that the physical neutral pion state is an admixture of the 
pseudoscalar octet fields 7r° and 77 ; that is, (7r°)phys = 7r° + ery with e ^ 0{md — rriu). Consequently e acts 
as an / = I “spurion” , encoding isospin-violating effects so that the matrix elements with the spurion 

are SU{2)f invariant. Isospin violation is also realized via the B+, mass difference; such effects are 
not encoded in the spurion framework, but they are also comparatively trivial. 

ii) Isospin violation can modify the form factors fi. The factor /o, e.g., is distinguished by the G-parity- 

violating decay to —> tt+tt”. The magnitude and phase of this effective p^-uj “mixing” can be elucidated 
from e+e“ ^ tt+tt” data however, the contribution is reflective of the decay B^ 7r+7r“7r°, 

so that ooo is modified in this region as well. In addition, electromagnetic effects distinguish y, probed 
in T decay, from /□ ||^, |^]. 

iii) Penguin contributions of |A/| = 3/2 character can occur, either through electroweak penguin effects [Q, 

or through isospin violation in the matrix elements of the gluonic penguin operator p^ . 

The impact of these isospin-violating effects can be redressed, at least in part. Eor example, in B ^ pir 
decay, the ^ 5 / 2,2 amplitude appears in the combination A^/ 2^2 + ^ 3 / 2,2 throughout Eqs. (|4|j|,^). Moreover, the 
two amplitudes share the same weak phase, to a good approximation. This emerges because, unlike K ^ mr 
decay no “|A/| = 1/2 rule” apparently exists in B ^ tttt decay — though B{B^ —> 7r+7r°) has yet to 
be conclusively determined p^. Recent theoretical estimates suggest that the magnitude of the ratio of the 
|A/| = 1/2 to |A/| = 3/2 amplitudes in B ^ tttt decay is roughly 0.3 p^, so that the |A/| = 5/2 amplitude is 
driven by an underlying |A/| =3/2 weak transition. Strong-interaction isospin violation acts in concert with 
a |A/| = 3/2 weak transition to generate a |A/| = 5/2 amplitude, whereas electromagnetism can generate a 
|A/| = 5/2 amplitude from a |A/| = 1/2 weak transition. The size of strong-interaction isospin violation is 
typified by the — rj mixing angle = V^{md — mu)/4(ms — rh) with m = {nid + m„)/2; we note that 
y ja ~ 1.45 p^, enhancing the extent to which and A^/ 2 p share the same weak phase. To the degree 

that this is true, the phenomenological parameters of Eq. include | A/| =5/2 effects as well. Thus we see 
that the single, crucial assumption of the isospin analysis is that the CKM factor accompanies |A/| = 1/2 

transitions exclusively, for in this case the weak phases of the ^ 5 / 2,2 and A^/ 2,2 amplitudes are identical. We 
have shown that isospin-violating contributions built on the |A/| =3/2 short-distance, weak transition do 
not impact the isospin analysis in B pn. However, non-|A/| = 1/2 penguin effects, be they electroweak 
penguin contributions or contributions consequent to isospin-violating effects in the hadronic matrix elements 
of |A/| = 1/2 operators, present a irreducible hadronic ambiguity from the viewpoint of this analysis. 

Empirical information on the all-neutral mode, ago, is essential to the extraction of a; however, it is possible to 
bound the strong-phase uncertainty using bounds on ooo and its CP-conjugate ago Under the assumptions 
we have articulated, the bounds on the hadronic uncertainty realized in i? —> p7r decay [ p^ are not modified by 
the presence of a |A/| =5/2 transition. 

An isospin analysis of B tttt decay also permits the extraction of sin(2a) from the mixing-induced CP 

asymmetry in i? —> tt+tt" In this case, in contrast, the |A/| = 5/2 amplitude cannot be combined with 
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the |AJ| = 3/2 amplitude. With we have 


5 +_ = A{B° 




) — -■^^1/2,0 + -^[^3/2,2 - ^5/2,2] , 


^00 = A{B'^ TT^TT^) — —-^^1/2,0 ~ Y^[^3/2,2 ~ ^5/2,2] 


6+0 = A{B^ 


TT^TT 


V3 

") = 


3/2,2 ' 


V2 


A^ 


5/2,2 


( 10 ) 


( 11 ) 


( 12 ) 


As the case of B —> p7r, three isospin amplitudes describe three empirical amplitudes, so that we expect our 
parametrization to persist in the presence of isospin violation. The lowest-dimension operators of the effective 
weak Hamiltonian generate transitions of |A/| = 1/2 and |A/| = 3/2 character, so that in the absence of isospin- 
violating effects in the hadronic matrix elements, Ai/ 2,0 ^ ^0 and A^/ 2,2 ^ 2 , and two amplitudes sufhce to 

describe the three transitions. The isospin analysis in H ^ tttt relies on the relation ( 6 -|_ — 600 )/V^— 6+0 — 0 [^ . 

The right-hand side of this relation is proportional to ^ 5 / 2,21 so that the required relation is broken by amplitudes 
of |A/| = 5/2 character. With such effects the isospin analysis can fail to determine the true value of sin(2a) 

The small B rate makes the full isospin analysis difficult to effect, so that bounds on the hadronic 

uncertainty in the extraction of sin(2a) have also been constructed |3^, ^ |3^. The presence of the 

|A/| =5/2 amplitude, as well as that of electroweak penguins, imply that the bounds can underestimate the 
size of the hadronic uncertainty |M. However, bounds which rely on the neutral B modes exclusively, such as 


Eq. (83) of Ref. |33 , contain the same linear combination of |A/| = 3/2 and |A/| = 5/2 amplitudes throughout 
— so that our arguments concerning the |AJ| = 5/2 amplitude in H > p 7 r decay are germane here as well. 
We conclude, to the extent the |A/| =3/2 and |AJ| =5/2 amplitudes share the same weak phase, that such 
bounds are insensitive to the |A/| = 5/2 amplitude and yield more reliable bounds on the hadronic uncertainty. 


3 Extension of the isospin analysis: Inclusion of the cttt channel 


The B —> (T 7 r channel has definite properties under CP, so that it can be included in the H > p 7 r analysis as 
well. Defining Oqq = A(H° ^ cnr^), we have 


- “00 


= T 


00 _ —ia 


P 


iOO 


(13) 


T™ and P™ are unrelated to the parameters of Eq. (||), so that we gain four additional hadronic parameters. 
However, more observables are present as well. Including the scalar channel, we now have = A{B^ —> 

7 r+ 7 r“ 7 r°) = /+ a.]_ + /_ a _ + fo aoo + fa- aoo> where fa- is the form factor describing a —> tt+tt”. It is worth 

noting that any discernable presence of the B —> uTr channel in the H > p 7 r phase space falsifies the notion that 
the “nonresonant” background can be characterized by a single, constant phase across the Dalitz plot For 
further discussion of the treatment of nonresonant contributions, specifically in ^ 37 r decay, see Ref. ]38| — 
note also Ref. 0. 

Neglecting the width difference of the B-meson mass eigenstates, as AT = T// — Tl and |Ar| ^ T = 
{rH + El)/ 2 , we note that the decay rate into 7 r+ 7 r“ 7 r° for a meson at time t = 0 is given by 


T{B^{t) TT+TT TT*^) = |A 3 ^pexp(-rt) 


l + |A3.n-P (1 ~ 


cos(Ato t) — Im sin( Am t) 


(14) 


whereas the decay rate into tt+tt 7 r° for a meson at time t = 0 is given 

1 + |A3.n-P (1 ~ l"^37rP) 


T{B^{t) ^ TT+TT TT*^) = |A 3 ^pexp(-rt) 


cos(Am t) -I- Im Aa^ sin(Am t) 


(15) 


We note that Asjr = ^Aa^/pAa^, where we have defined Aa^ = A(P° ^ 7 r+ 7 r“ 7 r°), and Am = Mh — Ml- 
Different observables are possible. For example, we can consider untagged observables, for which the identity 
of the B meson at t = 0 is unimportant, so that r(i3°(t) ^ 7r+7r“7r°) -|- r(i?°(t) ^ 7r+7r“7r°) oc (1 -f jAaTrP), or 
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we can consider time-integrated, tagged observables, containing 7r+7r“7r°) — 7r+7r“7r°), 

which are sensitive to (1 — |A 3 ,n-p). The products fif* contained therein are distinguishable through the Dalitz 
plot of this decay and thus the coefficients of these functions are distinct observables |^] . Were we to neglect the 
an channel, nine distinct, untagged observables exist, so that all the hadronic parameters save one would be 
determinable from the untagged data, for which greater statistics will be available | W . If we enlarge the analysis 
to include the an^ channel, the additional interferences possible imply that there are now sixteen distinct, 
untagged observables. Moreover, there are fifteen, rather than eight, tagged, time-integrated observables as 
well. Nevertheless, it would seem that the additional hadronic parameters associated with the cr7r° final state 
can be extracted from untagged data alone. Of course the practicability of the procedure relies on the amount 
of data eventually collected; moreover, the observables are highly correlated. 


4 Evaluating B —^ tt^ decay 

The effective, |Ai3| = 1 Hamiltonian for b dqq decay is given by 


HeS 


V2 


10 


A„(CiO“ + C20^) + Ae(CiOJ + C20^) - At y ao. 


i=3 


(16) 


where Xq = with an element of the CKM matrix. The Wilson coefficients Ci and operators Oi are 

detailed in Ref. Q, though we shall interchange CiO^ ^ C 2 O 2 so that Ci ~ 0(1) and Ci > C 2 ^ C' 3 ...io- 
The contributions with i = 3... 6 correspond to strong penguin effects, whereas those with i = 7... 10 are 
characterized by electroweak penguin effects. The Wilson coefficients with i = 7... 10 are numerically smaller 
than those with i = 3... 6, and penguin effects are not CKM-enhanced in 6 ^ dqq decay, so that we shall 
neglect the terms with i = 7... 10 all together. 

The decay amplitude for B M 1 M 2 , where Mi and M 2 are mesons, is given by 


A{B M1M2) = {MiM2\Hes\B) . 

The requisite matrix element contains terms of the form 


(17) 


C,{ii){MiM2\0,\B) . 


(18) 


We adopt the naive factorization approximation to effect estimates of the hadronic matrix elements. To wit, 
we separate Oi into a product of factorized currents, ji ® j 2 , and evaluate (Mi|ji|il)(M 2 |j 2 | 0 ), so that the 
operator matrix element becomes a product of a form factor and a decay constant. Such a treatment, albeit 
simple, is incomplete. The amplitude A{B M 1 M 2 ) is related to a physical observable and as such must be 
/x-independent, though the Ci therein do depend on /x. Evidently the /x dependence of the operator matrix 
elements compensates to yield a /x independent result. In the naive factorization approximation, we have 
replaced the operator matrix element by a product of a form factor and decay constant. These quantities are 
themselves physical observables and thus are without /x dependence, so that the overall /x dependence of the 
computed amplitude remains. Effecting this approximation, however, allows us to realize a clear connection 
to earlier work |^, |^, for our purpose is to illustrate the impact of using a scalar form factor consistent with 
low-energy constraints. 

In the naive factorization approximation, we can replace the effective Hamiltonian by the sum of products 
of factorized currents, so that HeS = + T®’®, where 


_ _ 

TA2 = -^A„ [ai?X7'"(l-75)&(8)d7^(l-75)u-ba2d7^(l-75)&(8'W7/x(l-75)w] , (19) 


T-3,4 _ \ 

^ “"Tf ‘ 


03 y d'y^{l - j 5 )b (g) 97^(1 - 75)9 + 04 y g 7 ''(l - 75)^ 0 dj^{l - -i^)q 


and 




Os y d - fF {\ - 75)6 g) 97^(1 + 75)g - 206 y ?(! - l 5 )b ® d(l + 75)9 


( 20 ) 


( 21 ) 
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We define Ui = Ci + Ci+i/3 for i odd and Ui = Ci + Ci-ij3 for i even; note, too, that q S u, d, s, c if ^ < rrif,. 

We now specifically consider B —> air transitions; the relevant formulae for B ^ pir are detailed in App. i By 
“tr,” we always mean a two-pion state with total isospin zero and in a relative S-wave state, (7r7r)s, understanding 
its dynamical origin in the strong pionic FSI for these quantum numbers — see the following section for a more 
detailed discussion. The matrix elements involving the a are 

q^{a{p,)\uj^il - 75)&|B°(pb)) = -*(M| - , (22) 


where q = ps — Pa and {a\q'^^{l + 75)q^|0), noting q' £ u,d,s, vanishes by C-invariance. It is just such a 
suppression mechanism that prompts the authors of Ref. |41 to argue that non-factorizable effects ought be 
relatively enhanced in two-body B decays to final states with scalar mesons, as the factorization contribution 
is itself small. Nevertheless, as our focus is the scalar form factor itself, we proceed with our estimates. 

For the ir'^ we have 

(7r+(p)|u7^(l - 75)d|0) = iUp>^ 

and 


(7r+(p)|M7^(l - 75)&|5 °(pb)) = 
With these relations, Eqs. 

(7r“cr|7deff|B~) = 


. , (M| - m2) 

{PB P)fi 2 


Ff—(g2) {Ml-Ml) 
2 




Grf 

V2\ 


and the equations of motion for the quark fields, we have 

aeMl 


/.(M2 - M2)F"--(M2) 


Autti ~ Ata4 -l- At 


m(mii + rh) 


+ At2a6 4^^^(M|-M2)Fo^--(M2) 


(mb - m) 


and 


{7r°a\n,s\B°) = 


Grif^ 


V2\V2 

— Xt2ae 


(M| - Ml)F^-%Ml) 

a\dd\0) 2 

' ' ' [M% - Ml) ° 


Xu(i2 + Ata4 ~ At 


06^2 


m{mi, + m) 


{mb — m) 


"(M2) ] 


V2 j 


(23) 

(24) 


(25) 


(26) 


where we have replaced the u, d quark masses with m and set = M^± = and M^i = Mgo = Ms, 
as we neglect isospin-violating effects. We adopt the usual phase conventions for the flavor wave functions 
7 r+, 7 r°, TT~ = ud, (uu — dd)/\/2, du, and adopt analogous relations for the p mesons as well. In the context of the 
R —> pTT analysis, the decay B^ an, specifically its penguin contributions, modify Eq. (^, the assumption 
crucial to the analysis. Thus it is important to make an assessment of the size of penguin effects in this 
decay. The terms containing uq, the scalar penguin contribution, are formally 1/mb suppressed, but can be 
chirally enhanced: the numerical factor /{mmb) 0.6 is only modestly less than unity. The second term 

proportional to qq contains (cr|(i(i| 0 ), where we anticipate, in the vicinity of the a resonance, ro. 7 r 7 r(s)(CT|d(i| 0 ) = 
( 7 r+(p+) 7 r"(p_)|Jd| 0 ) = y^TSRoFi *(s), where Fj‘(s) is the scalar form factor, which we detail in the next 
section, and s = {p+ +p-)'^. The parameter Bq is related to the vacuum quark condensate. Neglecting small 
terms of second order in the quark masses, Bq = — {0\qq\0)/F^, where F^-, the 7 r° decay constant, is 
Commensurately, we can simply set Bq = M^/(2m) to realize our numerical estimates. Note that describes 
the a —> 7 r+ 7 r“ form factor. With our conventions, Bq > 0, so that the two uq contributions are of the same 
sign. Using the parameters of Ref. and F^ttii- = r”*x with, as we shall determine, y = 20.0GeV“^, we 
find that the oe term containing to be roughly a factor of four larger. The {a\dd\0) term, present in the 

penguin contributions in B an, slightly enhances the qq contribution and its subsequent cancellation with 
the 04 contribution, as 04 and oe are of the same sign. The same cancellation occurs in B^ —> p^n^ (see App. 
A for a compilation of the relevant formulae), so that penguin contributions in > an^ can be expected to 
be crudely comparable to those in F° ^ p° 7 r°. Modifications of Eq. (^ can thus be expected to occur. 

In the previous section we determined how the isospin analysis in F —> p 7 r could be extended to include 
B an. Eqs. (|2|) and (^^, however, are related by isospin symmetry, so that it is useful to determine whether 
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additional constraints on in Eq. ( |l^ can be realized. Parametrizing the amplitudes in terms of 
we have 

A{B^ cr7r°) = - -^^3/2,1 (27) 

and 

A{B^ CTTr’*") = —Ai/2,i — 2^3/2,! , (28) 

so that two isospin amplitudes appear for the two empirical amplitudes. Although the |AJ| = 1/2 amplitudes 
are simply related, we see that no useful constraint emerges, as the tree amplitude in B^ —> (T7r+, which includes 
|A/| = 1/2 and 3/2 amplitudes, gives rise to two additional, undetermined hadronic parameters. 

We have illustrated that the penguin relation essential to the isospin analysis in i? —> pir, Eq. (||), can be 
broken through the presence of the B —>■ air decay channel. In our numerical estimates, however, we neglect 
penguin contributions, in order to retain a crisp comparison with earlier work [^, for our purpose is to illuminate 
the impact of the a —>■ tt+tt” form factor. With such an approximation, a computation of the Wilson coefficients 
in leading order in ag suffices Q, so that Ci{fj,) = 1.124 and C 2 {p) = —0.271 at fi = mt = 4.8 GeV as per 
Ref. [Q, to yield oi = 1.034 and 02 = 0.104. In contrast. Ref. ^ uses the “fitted” values C'i(mf,) = 1.105 
and C 2 {mb) = —0.228, to yield oi = 1.029 and 02 = 0.140. These are very similar to the Wilson coefficients 
in next-to-leading-order QCD, after Ref. Q, used in the B ^ pir analysis of Ref. Q. Eor definiteness, we 
shall adopt these last values, detailed in Sec. in our numerical analysis. The values of oi are quite similar, 
whereas those of 02 differ by tens of percent. Generally, we expect our numerical predictions for decay channels 
controlled by 02 to be less robust, as the scale dependence of 02 , as illustrated in Table 111 of Ref. Q, is severe. 
Note that it persists to a significant degree in the next-to-leading order treatment of Ref. as well. 

Nevertheless, let us proceed to consider numerical estimates for B —> air decay. We reconstruct the a meson 
from the (7r+7r“)s final state, so that we have 

Acr{B 7r+7r“7r) = ((cr 7r+7r“)7r|7Yeff|S) = A{B aTT)Tc„vnr , (29) 


where is the a —> tt+tt” form factor introduced previously. The scalar form factor contains the meson loop 
function and thus a regularization scale pr at which it is evaluated. In our approach, this scale dependence is 
disjoint from that associated with the renormalization of the operators of the effective, weak Hamiltonian, so 
that it is chosen for convenience and is quite independent of p. It may seem untoward to graft two very different 
calculations, namely, of A{B —s- utt) and of To-TrTr, to yield A{B —> (ct ^ 7r+7r“)7r). In a holistic treatment 
one might hope to recast a resonance and its subsequent decay products in terms of a single, complex hadron 
distribution amplitude, 4){x,p), which describes the non-perturbative dynamics. The analysis of the decay 
amplitude could then proceed via standard pQCD techniques In this manner the pr dependence of 

which we have spoken is connected, albeit loosely, to the scale dependence of cj){x,p). Nevertheless, the explicit 
“QGD factorization” analysis of Ref. shows that the scale dependence of the hadron distribution functions 
is trivial in an NLO analysis in ag , so that the consistency issues to which we have alluded are beyond the scope 
of current calculations. 

Turning to B~{pb) 7r+(p+)7r“(pi)7r“(p2) decay, we define u = (j)b — PiY = {P+ + 7 * 2 )^ and t = (p+ -|- 
PiY■ The contributions driven by the a resonance are of the form B~(pb) (cr ^ 7r+(p+)7r“(pi))7r“(p2) or 
B~(pb) ^ (cr —> 7r+(p+)7r“(p2))7r“(pi) — the latter is illustrated in Fig. y. (The corresponding contributions 
to B ^ pTT decay are illustrated in Fig. ||.) The two contributions add coherently, so that the branching ratio 
for B~(pb) —r 7r+(p+)7r“(pi)7r“(p2) is enhanced through the presence of two identical pions in the final state. 
In B° decay this does not occur, and we have B^{pb) {p+)tt~{pi)tt^{ p 2 ). Thus we find 


((^ ^ ^+^-) I Heff I B-) = a^ [r...(t) + , 


iB^<y)fnx2\/Ax2 A^2^ A 




In Ref. 1^, the cr ^ tt+tt vertex function is chosen to be 

^ cnriri^') — da'K+'n 


X — + iT^(x)M^ 


(30) 

(31) 


(32) 



where the running width ro-(x) is defined as 


r„(x) = 




(33) 


We shall adopt, however, the definition 

=xrr(a^) , (34) 

where the normalization x is fixed to be identical to that of Eq. (^), namely 


X 




(35) 


Using the parameters , Ma, and Ecr of Ref. ||], that is, Ma = 478 ± 24 MeV and Ecr = 324 ± 41 MeV as 

reported by the E791 collaboration in £)+ ^ 37 r, and (/o-TTTr = 2.52 GeV, we find x = 20.0 GeV“^. Alternatively, 
X = \/2/3ilo/(o’|d(i|0), so that this procedure determines (cr|Jci|0) as well. 

The cr meson is a broad I = J = 0 enhancement, close to the p meson in mass, so that B —> air decay can 
contribute to the allowed phase space of B ^ pir decay as well. To ascertain the impact of the B ^ an channel 
to B pn decay, we combine the decay channels at the amplitude level and then integrate over the relevant 
three-body phase space to determine the effective B ^ pn branching ratio. Eor B{pb) {p+)n~ {pi)n(p 2 ) 

decay, we define 


cos 6 — 


P'l ■ P 2 

|p'll|P2l ’ 


(36) 


where the primed variables refer to the momenta in the rest frame of the 7r+ {p'ff}n {p'l) pair, so that -l-p(|_ = 0. 
Letting P 2 = M|, we have 


E (R ^ PeffTr) = f dcosO f 

J-i J(j 




dt 


(Mp-sy 32(47rMB)3 




(37) 


where “peff” is determined by [Mp — Sff < t < {Mp + /3' refers to the velocity of particle i in the primed 

frame, and A4 is the sum of the amplitudes of interesij^. Typically <5 ~ (1 — 2)Ep. 

The scalar form factor we adopt describes the /o(980) ^ tt+tt” vertex function as well, albeit with a 
new normalization factor, as fa^ 7 ,{t){fo\dd\Q) = {n+{p+)TT (pi)|dd|0) = y/^BoT'^*{t), in the vicinity of 
the /o(980) resonance. The appearance of the /o(980) resonance in the tt+tt” channel is complicated by the 
opening of the KK threshold; the /o(980) ^ tt+tt” form factor is decidedly not of Breit-Wigner form. It should 
be noted that the unitarization procedure we employ here neglects the pp channel, though, as we discuss in 
the next section, this has no impact on the description of tttt scattering below s ~ 1.1 GeV. Note, too, that 
although multiparticle final states, particularly the 47r state, can contribute, they are demonstrably small for 
s < 1.4 GeV Q. The form factor(s) we adopt can be tested though the shape of the /o(980) contribution 
in B ^ /o(980)7r —> 37r, as well as through that in B ^ /o(980)7r — s- 


5 Coupled-channel pion and kaon scalar form factors 

In the preceding section, we encountered the non-strange, scalar form factor of the pion, E"(t). Such scalar 
form factors play a unique role in strong-interaction physics because they measure the strength of the quark 
mass term = muuu + ruddd -|- ..., i.e., the explicit chiral-symmetry-breaking term in QCD. However, 

since no scalar-isoscalar sources exist, these form factors cannot be determined directly but, rather, must be 
inferred indirectly from hadron-hadron scattering data. The most prominent example in this context is the pion- 
nucleon sigma term, which can be extracted from the analytically-continued and Born-term-subtracted isoscalar 
S-wave amplitude. Similarly, the pion scalar form factor can be obtained from meson-meson scattering data, 
with the added complication of channel couplings above the KK threshold at ffs ~ 1 GeV. At lower energies, 
the pion scalar form factor can be calculated in CHPT at one- Q and two-loop accuracy MB- The one- 
loop representation fails at surprisingly low energies, a consequence of the strong pionic ESI present in this 

^®For completeness, we note that u = ^)i with = (t + ~ Xfj )/(2\/7), E[ = (t + Mf — 

Xt+)/(2\/7), and Ej = (tttg — t — M|)/(2\/t), where and pj = Mf. 


9 









channel |^. This is also signalled by the very large pion scalar radius, (r|)^ ~ 0.6fm^, which is sizeably bigger 
than the pion vector radius governed by the rho mass ~ 6/M^ ~ 0.4fm^, indicating a smaller breakdown 

scale in the scalar-isoscalar channel. In fact, these FSI are so strong that they can dynamically generate the a 
meson, with no need of a genuine (pre-existing) quark model state . Furthermore, extending such an analysis 
to three flavors, one finds that all the scalar mesons with mass below 1.1 GeV can be dynamically generated, 
so that the genuine quark model nonet would have its center of gravity at about 1.3 GeV, see e.g. Ref. Q 
and references therein. Irrespective of this admittedly controversial assignment of the scalars, for further recent 
discussion see Ref. |^, the pion scalar form factor cannot be represented simply in terms of a scalar meson 
with a certain mass and width. Although the description of the low-lying scalar states in terms of dynamically 
generated, rather than “pre-existing,” states may be controversial and thus subject to ongoing discussion, let 
us stress that the scalar form factor itself is not. The form factor which emerges for the chiral unitary approach 
adopted here is quite comparable to that which emerges from the dispersion analysis of Ref. . Gonsequently, 
when one has a source with vacuum quantum numbers coupled to a two-pion state, it can be very misleading 
to use a simple Breit-Wigner parametrization, albeit with a running width. We shall display this graphically 
after discussing the scalar form factor and its construction, to be used in calculating the branching ratios to 
which we have alluded. We note in passing that this casts doubt on the extraction of the a meson properties 
from D ^ 3tt decays, note Ref. j^. This is quite in contrast to the pion vector form factor and the p meson, 
for which such a description works to good accuracy. 

We now summarize how to calculate the scalar form factor in Eq. (Q), following Ref. 0 , to which we refer 
in all details. Of course, one could also use the numerical result of the dispersion analysis of |0| for the scalar 
form factor. Note that the band found there overlaps with the band of Ref. ||l^ if one allows for parameter 
variations within known bounds. As a first step, we prefer to work in the framework of the chiral unitary 
approach, as it yields a form factor which is more convenient for applications. As we have noted, we need the 
scalar form factor of the pion for momentum transfers larger than 1 GeV. Consequently, we cannot work with 
the scalar form factor as computed in CHPT, but must invoke some resummation technique, as well as account 
for the channel coupling between the tttt and the KK systems. The resummation method is constrained only by 
unitarity and thus is not entirely model-independent. However, it can be strongly constrained by requiring that 
the so-constructed form factors match to the CHPT expressions, in the region where CHPT is applicable. Due 
to the channel coupling, we have to consider transition matrix elements of the non-strange and strange scalar 
quark bilinears between two meson states of isospin zero, namely tttt and KK, and the vacuum. We exclude the 
rjr] channel, as it does not affect the phase shifts or pion/kaon decays and transitions — it only plays a role in 
describing the inelasticity of / = J = 0 tttt scattering above 1.1 GeV, the physical rjrj threshold. For a detailed 
discussion of this point, we refer to Refs. |^, 57, The pertinent matrix elements are given in terms of four 
scalar form factor^^, 

(0|hn|7r7r) = V2BoT^is), {0\nn\KK) = V2 BoP^is) , 

0\ss\KK) = V2BoTf^{s) , 


SS TTTr) 


= V2Hor?(s) 


(38) 

where Bq is a measure of the vacuum quark condensate, Bq = — {0\qq\0)/, with iV ~ 93MeV the pion decay 
constant (strictly speaking its value in the limit of vanishing quark masses). In Eq. (^^ , the following notation 
is employed. The superscript s(n) refers to the strange/non-strange quark operator whereas the subscript 1(2) 
denotes pions and kaons, respectively. Furthermore, fin = {uu -\- dd)/^/^. The pion scalar form factors r”(s) 
and Ff (s) are calculated in Ref.[^ ^ through one loop in GHPT and r"(s) through two loops in Refs. |^, p^ . 
The scalar kaon form factors at next-to-leading order in GHPT were first explicitly given in jl^. Since it is 
central to our discussion, we give the explicit expression for r"(s) 


rr(s) = 


l + + ^ - LI) + 8(2Lg - 


3Ml 


F 2 


F 2 


+ /(s) + o /('5) f ’ 


with /(s) and /(s) given by 

^ 2s-Ml - ^ ^ s - , ^ M2 - ^ ^ 4s 

/(s) — 2^2 '^7r7i-(s) 0^2+ ^2 


2567r2 


M 


Ml 


4 log ^ - log —-f 3 




7 ( \ ^ ^ 1 ^ 7 ^ 12s / 1 

/(s) - 4 p 2 -^KKis) + T^^(s) + I L4 250^2 


log# + l 


(40) 


#®Here “tttt” and “XA'” states denote the linear combinations of physical tttt and KK states, respectively, with zero total isospin. 
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Here, Jpp(s) (P = TT,K,r]) is the standard meson loop function |^, and in this section, is the scale of 
dimensional regularization. The quantities /ip in Eq.(|^) are given by 




TOp 

327r2p2 



(41) 


Furthermore, the are scale-dependent, renormalized low-energy constants. We use here the same values 

as in Ref. [^. Since the combination 2Pg — L\ multiplies the normalization of r"(s) is sensitive to 

the precise value of these low-energy constants, which are only poorly known. This motivates our choice 
for the normalization of To-Trir)®)) given in Eq. (|35|). The one-loop representation can not be trusted beyond 
■y/s ~ 400 MeV, as it begins to diverge from the form factor extracted from tttt scattering data. To go to higher 
energies, one therefore has to study the constraints that unitarity imposes on the scalar form factors. The 
imaginary part of any scalar form factor is given by the appropriate meson-meson scattering T-matrix, so that 
the starting point for any unitary resummation scheme is exactly this scattering T-matrix |6Cf| , 

T{s) = [I + K{s)-g{s)]-^ -Ki^s) , (42) 


where s denotes the center-of-mass energy squared and K{s) can be obtained from the lowest order CHPT 
Lagrangian, e.g., K{s)ii = (s — M'^l2)lF'^. This T-matrix not only describes meson-meson scattering data but 
also, after gauging, photon decays and transitions, as reviewed in Ref. ||I^. In Eq.(|4^), the diagonal matrix 
g{s) is nothing but the familiar scalar loop integral 


9{s)t 


1 

(47r)2 


-1-f log^ 

/i2 


7'i(s) log 


g,(s) -I- 1 \ 

(Ti(s) - 1/ ’ 


(43) 


given here in dimensional regularization for the MS scheme, and cri(s) = -^/l — AMf /s. In what follows, we 
will set the regularization scale g, = 1.08 GeV. Of course, observables do not depend on this choice, and we 
could choose another value for g. However, the original investigation of meson-meson scattering transitions 
and decays by Oiler and Oset uses a three-momentum cut-off in the pion loop function. Translated to 
dimensional regularization, this gives the stated value of g. For energies above the threshold of the state i, 
unitarity implies the following relation between the form factors and the isospin-zero scattering T-matrix, 


Im r(s) =T(s) • .r(s) , 


STT-yi 


employing an obvious matrix notation with 


0 P2(s)0(s-4M|) 


r(s) 



(44) 


(45) 


where Pi{s) = \/sjA — Mf is the modulus of the c.m. three-momentum of the state i. Substituting Im r(s) 
with (r(s) — r(s)*)/(2z) and T(s) with the expression of Eq. ( ^^ and using the properties of the matrices g{s) 
and K{s), one can express r(s) as: 

T{s) = [lFK{s)-g{s)]-^ -R^s) , (46) 


where R{s) is a vector of functions free of any singularity. We remark that this procedure of taking the final 
state interactions into account is based on the work of Ref. |^. We also wish to stress that Eq.(^) can be 
applied to any K-matrix without unphysical cut contributions. As the final step, one fixes the functions in 
the vector R{s) by requiring matching of Eq.(^^ to the next-to-leading order (one loop) CHPT tttt and KK 
scalar form factors. This matching ensures that for energies where CHPT is applicable, these form factors fulfill 
all requirements given by chiral symmetry and the underlying power counting. The matching procedure thus 
determines the vector i?(s), as detailed in Ref. j^. For completeness, we give the expression for i?" pertinent 
to the scalar form factor P" (s). 


R^{s)i 


/Sf 4(T£ + ^ 

V2|^+ T2 


: -b 


8(2Tg - LI) 

tp2 


■(2MI + 3M1) 


Ml 

32 7r2 T2 



, (47) 


II 













where we have used the Gell-Mann-Okubo relation 3M^ = 4M|. — M^. This representation of the scalar form 
factors is valid from threshold up to energies of about 1.2 GeV. This range could be extended to higher energies 
by including multi-particle states. In fact, as we have noted before, in the u-channel one encounters larger 
values of ^/s. Therefore, we simply match our representation at = 1.2 GeV to the following asymptotic 
forms 

Re r?(s) ^ , Im r?(s) ^ , s s- oo . (48) 

s 

We have checked that the final results are insensitive to this choice of the matching point. The asymptotic form 
of the real part of the scalar form factor follows from quark counting rules in the crossed channel; it has 
also been found in the dispersion analysis of the Higgs decay into two pions |56|| . For continued contact with 
Ref. 1^ we match to the asymptotic form of the vertex function of Eq. a more precise treatment, if it were 
warranted, would involve smoothly letting K{s) —> 0 and solving for the form factors in a manner consistent 
with unitarity — for further discussion, see Ref. [Q. In Fig. ^ we display the scalar form factor r”(s) for 
yi < 1.2 GeV. In the low-energy region, the modulus of the form factor has it maximum at ^/s ~ 0.46 GeV, 
very close to the central value of the a meson mass deduced by the E791 collaboration from analyzing the 
I? ^ Stt data. In fact, since we wish to examine the consequences of using the more general vertex, Eq.(^J), 
as compared to the choice Eq.(^), used by Ref. |^, we have hxed the normalization constant x such that 
Re r"(s) has the same value as the Breit-Wigner representation at y/s = 0.478 GeV. This amounts to setting 
X = 20.0 GeV“^. The peak at y/s ~ 1 GeV is due to the /o(980) and the opening of the KK channel. Also shown 
in Fig. 0 is the scalar form factor generated using the Breit-Wigner form with a running width, adopted in 
Refs. H,^. The differences between this form factor and that deduced from the low-energy effective field theory 
of QGD are marked. In particular, in Fig. ^ we see that the Breit-Wigner representation of Ta-wn is deficient in 
that i) ImrCT 7 r 7 r(s) has a different shape as s approaches physical threshold, s —> 4M^ and ii) ReFcrTr,!-)®) does 
not possess a unitarity cusp at s = It is also of the wrong sign in this limit. Moreover, the shapes of 

the two form factors are very different above y/s ~ 0.5 GeV — this has particular consequence for the B ^ pir 
analysis. 


6 Vector form factor 

Thus far we have considered the a meson contribution to B decay. In this section we turn to the p 

meson contribution, B{pb) ^ pn ^ 7r+(p+)7r“(pi)7r(p2), which presumably dominates for t ~ In analogy 
to Eq. (^, the amplitude for B decay as mediated by the p resonance, Ap{B —> is the 

product of a R —> pTT amplitude and a p —*■ tttt vertex function, We give the relevant formulae for 

B ^ pn ^ 7r+7r“7r decay in App. ^ In this section we focus on the construction of Fp^rir, generating a form 
which is consistent with all known theoretical constraints. We detail our procedure, as the form of Tp-^j^ is 
important to the goals of the B pir analysis: it drives the size of the interference between p states produced 
in different regions of the Dalitz plot. Our vertex function differs from that of Ref. jl^, as the latter adopt a 
Breit-Wigner form. We also compare our form with that adopted in Ref. |^]. At this point we should mention 
that the numerical differences are not large — simply because the pion vector form factor can be described fairly 
well by a Breit-Wigner form. We also note that a unitarized version of the vector form factor starting from 
tree-level GHPT, including resonance fields, has been presented in Ref. [^, based on the methods described in 
Sec. ||. It could be used equally well as the parametrization employed here. 

The vector form factor Fp{s) can be directly determined from e+e“ —> tt+tt" data in the p resonance region. 
General theoretical constraints guide its construction: charge conservation requires Fp{0) to be unity, and time- 
reversal invariance and unitarity lead to the identihcation of the phase of Fp{s) with the Z = 1, / = 1 tttt 
phase shift, (5j(s), in the region where tttt scattering is elastic, s < Moreover, Fp{s) is an analytic 

function in the complex s plane, with a branch cut along the real axis beginning at the physical threshold 
s = AM^. Below the two-pion cut at s = 4Mj, the vector form factor is real. Furthermore, at small s the 
form factor can be computed in GHPT, as detailed in Refs. |^, All these constraints are captured by the 
Muskhelishvili-Omnes (MO) integral equation |Q. For s < its solution can be written @0 

_ F,(s) = P{s)a(s) , (49) 

^"^The solution of the MO equation with inelastic unitarity, important for s > has been discussed in Ref. ^3. 
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where P{s) is a real polynomial and the Omnes function, n(s), contains all the phase information, 


n{s) 


tan (s) 


exp 


s da 01 (sO 

TT s' s' - a - ie 


ImFp(s) 

ReFp{s) 


tani5i(s), 


(50) 


where i5((s) is the phase shift of / = 1, L = 1 scattering. In the Heyn-Lang parametrization |^, 0(s) is 
approximated by the quotient of two analytic functions, which contain polynomial pieces and the one-pion-loop 
expression for the p self-energy function. P{s) is chosen to be of third order in s in Ref. |^. We use here a 
recent update of the pion form facto r based on the Heyn-Lang parametrization. Specihcally, we use the 
parameter set of “solution B” of Ref. [|21|, reflecting a fit to the e+e” ^ tt+tt" data in the elastic region, subject 
to the constraint that the model reproduces the empirical tttt scattering length in the / = 1, L = 1 channel, 
a) = (0.038 ± 0.0.002) ||^. In what follows, we neglect the presence of the u resonance, or effectively 

— uj mixing. The latter is an important isospin-violating effect visible in the e+e“ ^ tt+tt” data in the close 
vicinity of s = — the fits of Ref. ||^ do include it. 

To realize the vertex function rp^,r(s), we define 


rp-TrTr)^) — 


-Fp{s) 
fPI 


where, as described in App. the electromagnetic coupling constant of the p meson, is fp^ 
0.001 GeV^ 1^. Ref. adopts a Breit-Wigner form for the vertex function, namely 


t^BW 

^ pTTTT 


(s) 


9p 

s - Mp2 + iTpMp ’ 


(51) 

0.122 ± 

(52) 


with the parameters Pp = 5.8 and Tp = 150 MeV — we use Mp = 769.3 MeV[^ . The two forms are compared in 
Fig. 1^ — the Breit-Wigner form offers a reasonable description of the vector form factor, though differences can 
be seen. In particular, the imaginary part of the Breit-Wigner form does not vanish below physical threshold, 
as it ought. This deficiency can be repaired by giving the Breit-Wigner form a running width, i.e.. 


pRW 

pTTTT 


(s) 


U{s) 


9p 


a-M^ + in(s) ’ 

f Pjs) 

V~s\piM^p)J ^ 


(53) 

(54) 


where p{s) = x/s/^ — M^. This form, modulo the proportionality constant, is adopted by Ref. 0. Fig. I 
compares Eqs. (|^ ) and ( p^ — the two forms are really very similar, though the forms differ slightly as 
s 4M^. The phase of the form factors, namely tan“^((Imrp^^)/(Rerp^^)), is plotted in Fig. Unitarity 
and time-reversal invariance dictates that the phase be that of / = 1, L = 1 scattering; the phase shifts from 
the data of Refs. @0 are shown for comparison. The forms of Eq. (|5T|) and Eq. (|^) confront the phase shift 
data nicely. The agreement of the latter form is a particular surprise, as it contains only two free parameters. 
Apparently the form of the imaginary part is precisely captured by the tttt branch cut, so that the phase is 
accurately determined by an arctan prescription to unitarize the amplitude. This is also realized in the CHPT 
analysis including an explicit p meson, for details see Ref. 

In the application to follow, we need to evaluate rp^^(s) for s > (M^r -|- , so that at ^/s = (M^r -|- M^) ~ 

920 MeV, the form factor of Eq. ( ^l| ) is matched to the Breit-Wigner form of Eq. ([^), yielding Re Fp{s) ~ a'/s 
and Im Fp[s) ~ for large s. We now turn to a discussion of our numerical results in p and tr-mediated 

i? —*■ 37r decay. 


7 Results and discussion 

Eirst, we must collect parameters. Eor crisp comparison with Refs. & we adopt the parameters used 
therein but reiterate them here for convenience. For meson masses and widths we use Mb = 5.279 GeV, 
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IVU = 139.57 MeV, Mp = 769.3 MeV, Tp = 150 MeV, M„ = 478 MeV, and = 324 MeV. We neglect the 
i?+, lifetime difference and use tb = 1.6 • 10“^^ sec. For quark masses, we use mj, = 4.6 GeV and rri = 6 
MeV. As for the CKM matrix elements, we adopt the Wolfenstein parametrization |Q, retaining terms of 0{X^) 
in the real part and of ©(A^) in the imaginary part, using A — 0.806, p = 0.05, p = 0.36, and A = 0.2196. 
For the Wilson coefficients, we use Ci = 1.100, C 2 = —0.226, Ca = 0.012, Ca = —0.029, = 0.009, and 


Cq = -0.033, after Ref. 


Ref. 


A, 


(B- 


. For form factors and coupling constants we use Fq {M^) = 0.46, after 
\m^) = 0.37, = 0.29, = ^2(92.4 MeV) ~ 131 MeV, and fp = 0.15 GeV^. 

Finally, we use Pp = 5.8 and = 2.52 GeV when using the form factors of Refs. §0- 

Let us begin by computing the branching ratios for R —> p7r and B o'k decay. Assuming two-body phase 
space, we use Eqs. (|A.2| ) and (A^) with Eqs. ( A.IC - O^ ), as well as Eq. (^. The results are tabulated in the 


first row of Table |1|. In the treatment of Ref. |^, the branching ratios of B ^ Mi M 2 decay and its charge 
conjugate are identical, even with penguin contributions, as no strong phase between the amplitudes of differing 
weak phase has been included. 


Table 1: Branching ratios (in units of 10“®) for B ^ pir and B ^ an decay, computed at tree level. 
The numbers in parentheses include penguin contributions as well, after Ref. |p^ . The hrst row of 
numbers compute the branching ratios using two-body phase space. Although Br(p ^ 7r“'"7r“) ~ 1, 
Br(cr ^ 7r+7r“) ~ 2/3, so that the numbers in brackets reflect the branching ratio times 2/3. The 
rows labelled “3-body” compute B ^ pn ^ Sn and B ^ an ^ n^n~n decay, integrating over the 
entire three-body phase space. “BW” denotes the use of the form factors of Refs |^, Eqs. (^, 
whereas “RW” denotes the use of the vector form factor of Ref. [Q, Eq. Finally, 

denotes the use of the form factors we have advocated. 



B^ ^ p 7r+ 

R® ^ P+TT 

RO ^ pOTT® 

B p^n 

B —>■ an 

B° an° 

2-body 

21.6 (21.0) 

5.96 (5.94) 

0.237 (0.308) 

4.74 (5.00) 

15.6 [10.4] 

0.147 [0.0982] 

3-body (BW) 

22.5 (22.1) 

— 

— 

4.11 (4.33) 

8.31 

0.0739 

3-body (RW) 

22.4 (22.0) 

— 

— 

4.08 (4.30) 

— 

— 

3-body (*) 

22.3 (21.9) 

— 

— 

4.03 (4.25) 

11.7 

0.108 


Proceeding to treat R —> p7r —> 37r and B —> an —> 7r+7r“7r decay, we foll ow t he p an d a intermediate states 
to their nn final states. We realize the transition amplitudes as per Eqs. (AM), (Ab), or (|^), and integrate 
over the three-body phase space as per Eq. (|3^), computing the integral in t over [2M,n., Mb — Mt^]. With this 
procedure, the branching ratios for B p~n~^, B p^n~, and B p^n^ become identical; we simply report 
the final result in the B p~n^ column. In treating B~ n^n~n~ decay, we divide the total rate by 1/2, 
to compensate for integrating over equivalent configurations. Noting that Br(cr ^ n'^n~) ~ 2/3, the quantity 
in brackets in the first row includes the factor of 2/3, for comparison with the three-body results. Gomparing 
the two-body branching ratios with those computed by integrating over the entire three-body phase space, it is 
evident that the branching ratios do not agree. The deviations can be attributed to both interference effects and 
finite-width effects. As an example of the former, both the diagrams illustrated in Figs. and ^ as well as those 
with Pi <-> p 2 , contribute to B~ n'^n~n~ decay. Glearly the interference of these diagrams is not included 
when the B~ p°n~ or B~ an~ process is treated as a two-body decay. As an illustration of latter, note 


that the couplings Pp and pa 
rates, namely 


are typically chosen so that they reproduce the p ^ nn and a 


3U = r(. 




) = 


1 


167rM^ 


[Ml - 


decay 


(55) 


and 


487rM? 


{Ml - ^Mlf/^\gp 


(56) 


For the meson masses and widths we have used, these formulae yield gcr-K'K = 2.53 and Pp = 6.03, respectively. 
Adopting these couplings in place of those used in Refs. |^, does reduce the discrepancy. Note that it is a 
“finite width” effect in that reducing the numerical width of the p or tr meson, in concert with Eqs. 
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reduces the discrepancy between the two- and three-body treatments. It is worth noting, however, for the 
physical values of the meson widths, that there is no one fixed coupling go-n-Tr or gp which removes the discrepancy 
entirely — the needed coupling in any given case is sensitive to the form factor chosen, as well as to the masses 
of the other particles in the final state. The former is apparent from a comparison with the vector form factor of 
Ref. 0, Eq. (|^), for which we use gp = 5.8 as well. The normalization issue of which we speak is particularly 
relevant for the comparison of theoretical branching ratios, for B VP decay, e.g., to experiment. It is 
present regardless of the form factor used. That is, in the case of the vector form factor we adopt, Eq. ( 13 ), the 
determination of fp-y can also be modified by finite width effects. The sign and size of the mismatch between the 
two- and three-body phase space calculations can be quite sensitive to the choice of form factor, as illustrated 
by the scalar case. We set the normalization of the form factor of Sec. denoted by in the Tables, to 
that of the form factor of Ref. |^, Eq. This is the only manner in which the parameters M^,ra- enter 

our analysis. Were we to determine the normalization so that the two- and three-body computations of the 
B° —!■ tTTr'^ —!■ 7r'''7r 7r° branching ratio yield identical results, the effective impact of the a in the pir phase space 
would be reduced by some 10%. 


Table 2: Effective branching ratios (in units of 10 for B ^ pn decay, computed at tree level. 
The numbers in parentheses include penguin contributions as well, after Ref. @]- The form factors 
are defined as in Table |3. 


S [MeV] (f.f.) 

B^ ^ p 7r+ 

B^ p+TT 

B° p°Tr° 

B p^TT 

TZ 

200 (BW) 

15.1 (14.7) 

4.21 (4.24) 

0.508 (0.497) 

3.50 (3.68) 

5.5 (5.1) 

300 (BW) 

16.4 (16.0) 

4.74 (4.76) 

0.918 (0.908) 

3.89 (4.10) 

5.4 (5.1) 

200 (RW) 

15.1 (14.8) 

4.19 (4.21) 

0.468 (0.463) 

3.49 (3.68) 

5.5 (5.2) 

300 (RW) 

16.4 (16.0) 

4.69 (4.70) 

0.835 (0.831) 

3.87 (4.07) 

5.5 (5.1) 

200 (*) 

15.3 (14.9) 

4.26 (4.28) 

0.473 (0.467) 

3.49 (3.68) 

5.6 (5.2) 

300 (*) 

16.4 (16.0) 

4.75 (4.76) 

0.865 (0.859) 

3.85 (4.06) 

5.5 (5.1) 

<5 [MeV] (f.f.) 

5° ^ p+TT- 

B'^ p“7r+ 

B° p°7r° 

B+ p°7r+ 

n 

200 (BW) 

15.1 (14.7) 

4.21 (4.15) 

0.508 (0.615) 

3.50 (3.68) 

5.5 (5.1) 

300 (BW) 

16.4 (16.0) 

4.74 (4.67) 

0.918 (1.02) 

3.89 (4.10) 

5.4 (5.0) 

200 (RW) 

15.1 (14.7) 

4.19 (4.13) 

0.468 (0.571) 

3.49 (3.68) 

5.5 (5.2) 

300 (RW) 

16.4 (15.9) 

4.69 (4.62) 

0.835 (0.935) 

3.87 (4.07) 

5.5 (5.0) 

200 (*) 

15.3 (14.8) 

4.26 (4.20) 

0.473 (0.576) 

3.49 (3.68) 

5.6 (5.2) 

300 (*) 

16.4 (15.9) 

4.75 (4.68) 

0.865 (0.963) 

3.85 (4.06) 

5.5 (5.1) 


In Table ^ we report the B —> pn ^ Sir branching ratios, computed in the manner of Ref. fl^j . Our numerical 
results differ slightly from theirs (note that B^ ^ B^ in their Table III). The upshot is that our estimate of TZ 
with penguin contributions is ~ 5.1, rather than the 5.5 they estimate. We show the branching ratios computed 
for differing vector form factors; these differing choices have little impact on the resulting branching ratios, or 
on TZ. 

We compute the B —> utt branching ratios in Table ^ In this case our computed branching ratios, for 
B~ —> tT7r“ decay, are a factor of two larger with the same form factors and parameters input, as our formula, 
Eq. (|30|), differs from theirs by a factor of \f2. Thus the impact of the a in the p°7r“ phase space is rather 
larger than that estimated in Ref. |^. Updating the scalar form factor to use what we feel is its best estimate, 
we find that the values of TZ are smaller still. Interestingly, the computed values of TZ are comparable to the 
empirical results, albeit the errors are large. [An additional contribution to the phenomenological value of TZ, 
realized through a diagram mediated by the meson, is proposed in Ref. [f^ .j 

Turning to R —> cttt*' decay, we see that the contribution of the a meson to B^{B^) —> pn decay is much 
smaller — with the scalar form factor we advocate, the effect is some 10%. Interestingly the a has a tremendous 
impact on B~ —> p^T:~ decay, and a relatively modest one on B^ —> p^i:^ decay. Let us emphasize that we have 
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realized our numerical analysis at tree level. It is the relative size of the penguin contributions in 
and 5° ^ /9°7r° decay which is of relevance to the isospin analysis to extract a. The presence of the air' 
state in the phase space can break the assumed relationship, Eq. (|^), between the penguin contributions in 
piT and thus mimic the effect of isospin violation — alternatively we can expand the pir analysis to include the 
(TTF channel. Nevertheless, we expect our estimates to be crudely indicative of the importance of these effects 
— quantitatively, however, differences may exist. It is worth noting that the and p^i^^ contributions can, 
to some measure, be distinguished. Certainly the and pi:^ contributions behave differently under the cut 
on the invariant mass of the tt+tt” pair, recalling Eq. (|^. Moreover, making a cut on the helicity angle 9, 
defined in Eq. (|^, ought also be helpful in separating the and a contributions. This is illustrated in Eigs. 
and |[ The p^ir contributions roughly follow a cos^(0) distribution, whereas the cftt contributions are quite flat, 
save for the bump resulting from the ro. 7 r 7 r('w) term in Eq. (^). Cutting on the helicity angle 9 should also 
help disentangle the contributions from some of the B* resonances, discussed in Ref. |16 . The contributions 
of B* resonances to the pir channels should be included in a more refined analysis, but they will not alter the 
conclusions drawn here. 


Table 3: Effective branching ratios (in units of 10 ®) for i? —s- ctt and B ^ pn decay, computed 
at tree level. The form factors are defined as in Table |^. 


5 [MeV] (f.f.) 

B —> UTT 

B ^ (p° -f cr)7r 

RO ^ ctttO 

R° ^ (p° -f cr)7r0 


200 (BW) 

2.97 

6.16 

0.0258 

0.516 

3.1 

300 (BW) 

5.17 

8.61 

0.0457 

0.940 

2.5 

200 (RW) 

2.97 

6.19 

0.0258 

0.475 

3.1 

300 (RW) 

5.17 

8.62 

0.0457 

0.855 

2.4 

200 (*) 

4.11 

7.61 

0.0396 

0.508 

2.6 

300 (*) 

7.01 

10.7 

0.0663 

0.916 

2.0 


8 Summary 

In this paper, we have scrutinized the role of the a meson in B pir 3 tt decay, understanding its dynamical 
origin in the strong pion-pion final state interactions in the scalar-isoscalar channel. The presence of the air^ 
contribution in the p^ir^ phase space is important in that it can break the assumed relationship between the 
penguin amplitudes, Eq. (|^, consequent to an assumption of isospin symmetry. In this, then, its presence 
mimics the effect of isospin violation. The salient results of our investigation can be summarized as follows: 

i) We have considered how SM isospin violation can impact the analysis to extract a in B pn decay. Under 
the assumption that |A/| =3/2 and |AJ| = h(2 amplitudes share the same weak phase, the presence of 
an additional amplitude of |A/| = 5/2 character, induced by isospin-violating effects, does not impact the 
B piT analysis in any way. This is in contradistinction to the isospin analysis in B tttt. Thus the 
isospin-violating effects of importance are those which can break the assumed relationship between the 
penguin contributions, Eq. (^. 

ii) The scalar form factor can be determined to good precision by combining the constraints of chiral sym¬ 
metry, analyticity, and unitarity. The form factor we adopt describes the appearance of the /o(980) as 
well, so that the shape of the /o(980) contribution in i? ^ /o(980)7r ^ 37r, e.g., should serve as a test 
of our approach. We emphasize that the resulting scalar form factor is very different from the commonly 
used Breit-Wigner form with a running width. This is in stark contrast to the vector form factor, which 
is dominated by the p resonance. In that case, one can construct simple forms that fit the theoretical and 
empirical constraints. 

iii) We have pointed out that the two- versus three-body treatments of the decays B — > pn, B ^ air can lead 
to differing results due to finite-width and interference effects. 
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iv) Remarkably, the impact of the cttf channel on the ratio TZ, cf. Eq. (H), is huge. The numbers we find for 
TZ are in agreement with the empirical ones, given its sizeable experimental uncertainty. This underscores 
the suggestion made, as well as improves the calculations done, in Ref. [||. Our analysis is based on 
consistent scalar and vector form factors. 

v) On the other hand, the impact of the air channel on the B ^ pir isospin analysis is merely significant. 
Varying the cuts on the tttt invariant mass and helicity angle 6 should be helpful in disentangling the 
various contributions. 

vi) We have shown that one can expand the isospin analysis to include the cttt channel because it has definite 
properties under CP. This may be necessary if varying the cuts in the tttt invariant mass and helicity angle 
9 are not sufficiently effective in suppressing the contribution from the air^ channel in the p^ir^ phase 
space. 

This work is merely a first step in exploiting constraints from chiral symmetry, analyticity, and unitarity in the 
description of hadronic B decays. In particular, the contribution of the “doubly” OZI-violating strange scalar 
form factor and its phenomenological role in factorization breaking ought be investigated. 
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A Formulae for B —^ pTT —^ tt^ decay 

In this appendix, we report the formulae needed to evaluate B ^ pn ^ transitions, as per the approach 

of Sec. For clarity of comparison, we conform as much as possible to the notation and conventions of Ref. [p^ , 
but give the formulae required for completeness. This also allows us to identify the changes in replacing the 
Breit-Wigner form adopted for the p resonance in Ref. [p^ with the pion vector form factor we discuss in Sec. 
Defining 

(7r°(p2)7r"(pi)|p"(pp,e)) = ppC ■ {jp2 - Pi) (A.l) 

and 


(p*(Pp)7r(p7r)|Weff|S°(pB)) = 2e*-p^r)\ 

{p°iPp)T^~{p7T)\'HeB\B~{pB)) = 2e*-p^Tf, 

where i G (+, 0, —), the B^ pn —> amplitude can be written 


(A.2) 

(A.3) 


Ap(B°(pb) ^ 7r+(p+)7r ipi)TT°{P 2 )) =-r]°(s - u)Tp^^{t) + (s - t)Tp^^{u) + rj {t - u)Tp^^{s) , (A.4) 

where the first of these contributions is illustrated in Fig. ^ We have used s = {pi +^ 2 )^, t = {p+ +pi)‘^, 
and u = (p+ +P 2 )^, and have summed over the polarization states of the p* mesons, setting M^± = Mt^o. 
With our conventions for the flavor content of the meson states, we see that |7r+), |p+) = —|/ = 1I^ = 1), 
whereas the other tt and p charge states do not have a minus sign when written in the isospin basis. Using the 
isospin-raising operator t+, we thus determine from Eq. (AU) that (p+)t^~( pi)\p^[Pp, e))_= —Pp e ■ {p+ ~Pi) 
and (7r+(p+)7r°(p2)|p^(Pp, e)) = 5p e • {p+ ~ P 2 )', the signs we indicate consequently follow^. For the B~ 
p^Tr~ —>■ TT'^TY~TT~ amplitude we have 


Ap{B (ps) ^ 7r+(p+)7r (pQtt {P2)) = -9° [{s - u)T p„^{t) + {s - t)rp^^{u)] . 


(A.5) 


Note that rp,r 7 r(s) is the pion vector form factor, for which a Breit-Wigner form is used in Ref. |0. As discussed 


Sec. we replace 


- pTTTT 


(x) = 


9p 


-Fp{x) 


X — + iTpMp 


(A.6) 


^®We thank J. Tandean for discussions on this point. 
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where fp^ is the electromagnetic coupling constant of the p meson, determined from 


r(p ^ e'''e ) 


47ra^ 2 


(A.7) 


where r(p ^ e+e ) is, in turn, extracted from e+e ^ tt+tt data at s = M^, as described in Ref. |Q. For 
the “solution B” ht of Ref. we have fp^ = 0.122 ± 0.001 GeV^ |^. The two forms are compared in Fig, [i. 
In Ref. 0 the parameters Pp = 5.8 and Fp = 150 MeV are chosen — we use the value Mp = 769.3 MeV[^ 
for the p meson mass, as it is not reported in Ref. |l^ . 

To determine rf, ff, we introduce 

(A.8) 

<l^{p~{Pp-,^)\uip{l-l^)b\B^{PB)) = -i‘2Mp{e* ■ , (A.9) 

where q = ps — Pp, and recall Eqs. (^ ) and (|^), to find 


+ _Gf 

V2 

Gf 


AuUl — XtO-4 + At 


aeM^ 


m{mb + to) 


UMpA^^P{Ml) , 


q — fpFi {Alp) 




2^/2 I . 


^u0^2 At(l4 — \t 


aaMl 


rh{mb + to) 


UMpA^-P{M'^) + [Xua2 + Xta4] fpF^^^iMl) 




XuCii ~ XtQi + At 


aeM^ 


m{mb + rh) 


UMpA^^P{M'^) + [Xua2 + Ata4] fpF^^^{Ml) 


(A.IO) 

(A.ll) 

(A.12) 

(A.13) 


where we neglect electroweak penguin contributions, as well as all isospin-violating effects. Our expressions 
agree with those of Ref. and Ref. . 
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FIGURES 



Figure 1: B ^ decay as mediated by the cr resonance. The factorized 

weak vertex is denoted by “(g)®”. The filled circle denotes the strong three- 
meson vertex, here cr ^ 27r. 



Figure 2: B ^ decay as mediated by the p resonance. The factorized 

weak vertex is denoted by “(g)®”. The filled circle denotes the strong three- 
meson vertex, here p —>■ 2 tt. 


21 




Figure 3: The a 7r+(p_|_)7r“ (p_) form factor Fo-TrTr as a function of yfs, with 
s = (p+ + p-Y. The real (solid line) and imaginary (dot-dashed line) parts of 
Fo-TrTr, as well as its modulus (dashed line), are shown. The curves which do 
not persist below physical threshold, ^/s = 2M,r ~ 0.27 GeV, correspond to 
the form factor adopted in Ref. |^, whereas the curves which extend to s = 0 
correspond to the form factor adopted here 
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Figure 4: The p {p+)'K~{p-) form factor —FpTrTi- as a function of \/s, 

with s = (p+ + p-Y. The form factor is shown in the region for which I = 
1, / = 1 TTTT scattering is elastic. The real (solid line) and imaginary (dot- 
dashed line) parts of Fp^^, as well as its modulus (dashed line), are shown. 
Noting Eq. (^), the arrows indicate the form factor given by Fp{s)/fp^, as 
detailed in Sec. whereas the other curves correspond to the Breit-Wigner 
form —gp/{s — + iMpTp), adopted in Ref. fl^ . 
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Figure 5: The p 7r+(p+)7r“(p_) form factor —Fp^.^. as a function of ^/s, 
with s = (p+ +p-)^. The form factor is shown in the region for which 1 = 1, 
/ = 1 tttt scattering is elastic. The real (solid line) and imaginary (dot-dashed 
line) parts of Fp,rir, as well as its modulus (dashed line), are shown. Noting 
Eq. (0), the arrows indicate the form factor given by Fp{s)/fp^, whereas the 
other curves correspond to the form of Eq. (^^ adopted in Ref. |^] . 
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Figure 6: The phase of the vector form factor rp,rir(s) as a function of ^/s, in 
the region where the scattering is elastic. The form factor we adopt, Fp{s)/fpj 
(solid line), the relativistic Breit-Wigner form of Ref. ||^ (dot-dashed line), as 
well as that of Ref. 0 (dashed line), are all shown. Unitarity and time-reversal 
invariance requires that the phase be the phase shift <5^ of / = 1, L = 1 tt-tt 
scattering. The empirical phase shifts of Ref. (□) and Ref. 0 (a) are 
indicated. 
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Figure 7: Absolut e squ are of the matrix element, |Mp, for B~ —> p°7r“ decay 
(dashed line), Eq. (AA), and for B~ a'iT~ decay (solid line), Eq. (^0|), as a 
function of cos 9 at t = M^. The scalar and vector form factors advocated in 
Secs. H and ^ have been used. The bump in the solid line reflects the presence 
of the rCT 7 r 7 r(w) term in Eq. 


0.8 


0.6 


0.4 


0.2 


0 



cos(0) 


Absolute square of the matrix element, |Mp, for B^ 

^0 


p°7r° decay 


Figure 8: 

(dashed line), Eq. (A^), and for ^ cr7r° decay (solid line), Eq. as a 
function of cos9 at f = M^. The scalar and vector form factors advocated in 
Secs. H and H have been used. 
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